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Abstract 


We present optical photometry and spectroscopy of SN 20161et (=Gaial6bvd=PS17brq), an unprecedented Type I 
supernova (SN I) at z — 0.0676 with no obvious analog in the existing literature. SN 2016iet exhibits a peculiar 
light curve, with two roughly equal brightness peaks (zz —19 mag) separated by about 100 days, and a subsequent 
slow decline by about 5 mag in 650 rest-frame days. The spectra are dominated by strong emission lines of calcium 
and oxygen, with a width of only 3400 km s_!, superposed on a strong blue continuum in the first year. There is no 
clear evidence for hydrogen or helium associated with the SN at any phase. The nebular spectra exhibit a ratio of 
Licam/Lio y = 4, much larger than for core-collapse SNe and Type I superluminous SNe. We model the light 
curves with several potential energy sources: radioactive decay, a central engine, and ejecta-circumstellar medium 
(CSM) interaction. Regardless of the model, the inferred progenitor mass near the end of its life (1.e., the CO core 
mass) is 255M, and potentially up to 120 Mc, clearly placing the event in the regime of pulsational pair 
instability supernovae (PPISNe) or pair instability supernovae (PISNe). The models of CSM interaction provide 
the most consistent explanation for the light curves and spectra, and require a CSM mass of =33 MG ejected in the 
final decade before explosion. We further find that SN 2016iet is located at an unusually large projected offset 
(16.5 kpc, 4.3 effective radii) from its low-metallicity dwarf host galaxy (Z ~ 0.1 Zo, L ~ 0.02 Ly, 
M x 10% Mo), supporting the interpretation of a PPISN/PISN explosion. In our final spectrum at a phase of 
about 770 rest-frame days we detect weak and narrow Ha emission at the location of the SN, corresponding to a 
star formation rate of 23 x 10 * Mo yr !, which is likely due to a dim underlying galaxy host or an H II region. 
Despite the overall consistency of the SN and its unusual environment with PPISNe and PISNe, we find that the 
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inferred properties of SN 2016iet challenge existing models of such events. 
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1. Introduction 


Over the past decade, wide-field optical surveys have been 
successful at uncovering new types of rare transient events, 
such as superluminous supernovae (SLSNe; e.g., Chomiuk 
et al. 2011; Quimby et al. 2011), tidal disruption events (e.g., 
Gezari et al. 2009; Arcavi et al. 2014; Chornock et al. 2014), 
fast blue optical transients (e.g., Drout et al. 2014), and 
calcium-rich transients (e.g., Kasliwal et al. 2012). At the same 
time, some predicted types of SNe remain elusive, in particular 
pair instability SNe (PISNe; Barkat et al. 1967) and pulsational 
pair instability SNe (PPISNe; Woosley et al. 2007). 

PISNe are the predicted explosions of massive stars with zero 
age main sequence (ZAMS) masses of 130-260 Mo. They are 
caused by the high core temperatures in massive stars, leading to 
the production of electron—positron pairs which in turn lead to 
contraction of the core, followed by explosive oxygen burning 
that eventually unbinds the star, leaving no remnant behind 
(Barkat et al. 1967; Rakavy & Shaviv 1967; Kasen et al. 2011). 
The luminous SN 2007b1 was claimed to be a PISN based on its 
slowly declining light curve and claimed large mass of 
radioactive “Ni (Gal-Yam et al. 2009) However, Nicholl 
et al. (2013) suggested that SN 2007bi was powered by a 
magnetar central engine with an ejecta mass of only 7 Mọ. The 


spectra, metallicity, and observational features of SN 2007bi 
have been shown to be inconsistent with a PISN (Dessart et al. 
2013; Jerkstrand et al. 2016; Mazzali et al. 2019). To date, there 
have been no definitive cases of PISN explosions. 

For stars with slightly lower ZAMS masses of 95-130 Mo, 
the progenitor can experience multiple non-destructive pair 
instability episodes that expel material prior to the final core 
collapse. These pulses can lead to shell collisions that power an 
SN-like transient—a PPISN—while the inner layers contract 
and temporarily resume stable burning (Woosley et al. 2007). 
The ejecta produced in the ultimate core collapse explosion (if 
successful) can also interact with shells ejected ~10* yr prior, 
powering a luminous and long-lasting SN (Woosley et al. 
2007; Aguilera-Dena et al. 2018). There are several observed 
SNe with unusual photometric or spectroscopic features that 
have been interpreted as signs of PPISNe. The hydrogen-poor 
SN 2010mb (Ben-Ami et al. 2014) exhibited a long-lasting 
light curve, narrow oxygen emission lines, and a persistent blue 
continuum. The hydrogen-rich iPTF14hls had a light curve 
with at least five distinct peaks over multiple years, absorption 
lines of constant velocity, and a sub-solar metallicity host 
(Arcavi et al. 2017). The hydrogen-poor SLSN iPTFl6eh 
showed evidence for a light echo from a shell whose properties 
might match some PPISN models (Lunnan et al. 2018). 
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Against this backdrop, we present detailed optical observa- 
tions of SN 20161et, a transient initially discovered by the Gaia 
Photometric Science Alerts Team, and also recovered by the 
Catalina Real-Time Transient Survey (CRTS) and the Pan- 
STARRS Survey for Transients. Our follow-up observations 
over the past three years place the event at a redshift of 
z = 0.0676 and reveal highly unusual light curves, spectra, and 
environment. The light curve has a double-peaked structure of 
roughly equal brightness peaks of ~ —19 mag separated by 
100 days, followed by a long duration and slow decline, and a 
total radiated energy of a few x 10% erg. The spectra are 
dominated by intermediate-velocity (=3400 kms) lines of 
calcium and oxygen superposed on an initially blue continuum, 
and show no clear evidence for hydrogen and helium 
associated with the SN. Finally, we find that SN 2016iet is 
located at an unusually large offset of 16.5 kpc (4.3 effective 
radii) from its low-metallicity dwarf host galaxy. Weak and 
narrow Ha emission at the SN position points to a potential 
underlying dwarf galaxy or a single HII region. The 
combination of these unusual properties prompts us to explore 
a wide range of models for the energy source powering 
SN 2016iet: radioactive decay, a central engine (magnetar or 
fallback accretion), and circumstellar interaction. We find that, 
regardless of the model, the inferred progenitor mass (CO core) 
is 255 M. (and perhaps as high as =120 Mo), placing this 
event in the realm of PPISNe and PISNe. 

The structure of the paper is as follows. In Section 2 we 
present the discovery of SN2016iet and our follow-up 
observations. In Section 3 we describe the multi-color and 
bolometric light curves, as well as the spectra. In Section 4 we 
present our modeling of the light curve, and in Section 5 we 
describe the host galaxy. In Sections 6 and 7 we discuss 
possible interpretations of this event and present our key 
conclusions. Throughout the paper we assume a flat ACDM 
cosmology with Ho = 69.3 kms! Mpc !, Q,, = 0.286, and 
Q4 = 0.712 (Hinshaw et al. 2013). 


2. Discovery and Observations 
2.1. Discovery 


SN 2016iet was first discovered by the Gaia Science Alerts 
(GSA; Wyrzykowski et al. 2012) on 2016 November 14 
(MJD — 57706.596) and announced two days later with the 
designation Gaialóbvd with a magnitude of G = 18.61 
at coordinates R.A. = 12°32"33523, decl. = +27%07'15"49 
(J2000). The transient was subsequently recovered by the 
CRTS (Drake et al. 2009) on 2017 January 5 with a magnitude of 
C — 19.03 and was given two separate designations (MLS170316: 
1232334-270716 and CSS170201:123233+270716), and by the 
Pan-STARRS Survey for Transients (PSST; Huber et al. 2015) on 
2017 March 6 with a magnitude of i = 18.72, and the designation 
PS17brq. The GSA and PSST provide upper limits from images 
prior to discovery of G > 20.7 on 2016 July 11 and i > 19.7 
on 2016 May 19; see Figure 1. We find a 3c upper limit of r 3 
23.0 mag through forced photometry at the position of SN 2016iet 
in archival Pan-STARRSI 37 images. 

As part of our search for exotic transients, we flagged 
SN 2016iet as a host-less transient given the PS1/37 non- 
detection at the SN position. We eventually found that it resides at 
a large separation of =11” from a galaxy with a half-light radius 
of =2"6 that we identify as its host based on a common redshift 
(Section 5). This prompted us to obtain an optical spectrum on 
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2017 April 30 with the Inamori-Magellan Areal Camera and 
Spectrograph (IMACS; Dressler et al. 2011) on the Magellan 
Baade 6.5 m Telescope. The spectrum revealed strong calcium 
and oxygen features, weak Si features, and no evidence for H or 
He lines, which in principle lead to a Type Ic designation. 
However, SN 2016iet does not resemble any typical SNe Ic, and 
we therefore adopt a general classification of an SN I. 

To determine the redshift of SN 2016iet we measure the 
wavelengths of the strongest emission lines in two individual 
spectra: [Car] AA7291,7324, Om 27320,7330, and Call 
AAA8498,8542,8662. We simultaneously fit a Gaussian profile 
for each emission line in each spectrum and derive z — 
0.0676 + 0.0002. 


2.2. Optical Photometry 


We obtained images of SN 2016iet in the gri filters from 
several instruments: KeplerCam on the 1.2 m telescope at Fred 
Lawrence Whipple Observatory (FLWO), Templeton on the 
1.3 m McGraw-Hill Telescope at MDM Observatory, the Low 
Dispersion Survey Spectrograph (LDSS3c; Stevenson et al. 
2016) on the Magellan Clay 6.5 m telescope at Las Campanas 
Observatory, and MMTCam on the MMT 6.5 m telescope. We 
processed the images using standard IRAF routines, and 
performed photometry with the daophot package. Instru- 
mental magnitudes were measured by modeling the point- 
spread function from each image using reference stars and 
subtracting the model from the target. For calibration, we 
estimated individual zero-points by measuring the magnitudes 
of field stars and comparing to photometric AB magnitudes 
from the PS1/37 catalog. The uncertainties reported here are 
the combination of the photometric uncertainty and the 
uncertainty in the zero-point determination. The resulting 
photometry is summarized in the Appendix. SN 2016iet was 
also observed on several epochs by Gaia, PSST, and CRTS in 
the G, i, and unfiltered C bands, respectively. These data were 
obtained from the public databases for the GSA System, the 
Pan-STARRS NEO survey, and the Catalina Sky Survey 
photometric data set, and are included in the Appendix. 

The full light curve of SN 2016iet, spanning to 2019 January, is 
shown in Figure 1, where we define phase O to be the date of the 
brightest magnitude in the G-band, MJD = 57732. We use this 
reference point throughout the paper, with all days quoted in 
the rest frame unless otherwise specified. All of the photometry 
has been corrected for Galactic extinction with Ay= 0.041 + 
0.001 mag (Schlafly & Finkbeiner 2011). We determine that host 
subtraction and host extinction corrections are not necessary, 
based on measurements of the host extinction being consistent 
with zero and the large separation of SN 2016iet from its host 
(Section 5). All the photometry collected for this work is made 
available on the Open Supernova Catalog” (Guillochon et al. 
2017). 


2.3. Optical Spectroscopy 


We obtained nine epochs of low-resolution optical spectrosc- 
opy at phases of 132-772 days. We used the IMACS and 
LDSS3c spectrographs on the Magellan telescopes, the Blue 
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Figure 1. Optical light curves of SN 2016iet in the gri, Gaia-G, and CRTS unfiltered-C bands. The photometry is corrected for Galactic extinction. The vertical dashed 
lines indicate the epochs of spectroscopic observations. The inset shows a zoomed-in on the initial double-peaked structure (marked by arrows), with the photometry 


binned for clarity. 


Channel (Schmidt et al. 1989) and Binospec (Fabricant et al. 
2003) spectrographs on the MMT, the Gemini-North Multi- 
Object Spectrograph (GMOS; Hook et al. 2004), and the Ohio 
State Multi-Object Spectrograph (OSMOS; Martini et al. 2011) 
on the 2.4 m Hiltner Telescope at MDM Observatory. All the 
spectra were obtained with the slit oriented along the parallactic 
angle. Details of the spectroscopic observations are provided in 
the Appendix. 

We processed the spectra using standard IRAF routines with 
the twodspec package. The spectra were bias and flat 
corrected, the sky background was modeled and subtracted 
from each image, and the one-dimensional spectra were 
optimally extracted, weighing by the inverse variance of the 
data. A wavelength calibration was applied using a HeNeAr 
lamp spectrum taken near the time of each science image. 
Relative flux calibration was applied to each spectrum using a 
standard star taken on the same night. The spectra were then 
calibrated to absolute flux using our FLWO gri photometry as 
a reference, where a normalization constant was applied to 
each spectrum to match the expected flux from the photometry 
using the PYPHOT Python package. Lastly, the spectra were 
corrected for Galactic extinction and transformed to the rest 
frame of SN 2016iet. All the spectroscopy collected for this 
work is available on the Open Supernova Catalog (see footnote 7); 
Guillochon et al. 2017). 

The spectra obtained on 2017 November 25 and 2018 
January 21 with the Blue Channel Spectrograph also include 
the host galaxy. The one-dimensional spectra of the host were 
extracted in the same way as outlined above and are discussed 
in Section 5. 


3. Observed Properties of the Light Curves and Spectra 
3.1. Multi-band Light Curve 


The multi-band light curves of SN 2016iet are shown in 
Figure 1. The most striking features are the two distinct peaks, 
and the overall long duration of the event. The first peak, 
detected by Gaia, appears to have a roughly symmetrical shape, 


rising by 0.5 mag in 30 days to a peak of Mg ~ —19.2, and 
declining to Mg ~ —18.7 after 25 days, when the SN was 
detected by the CRTS and PSST. These data reveal a second 
peak, with Mc = —18.7 mag, about 100 days after the first peak. 
We stress that the distinct peaks are not an artifact of comparing 
different filters, since the rise and decline of the first peak is 
observed entirely by Gaia, while the rise and decline of the 
second peak is clear even from the CRTS data alone. 

Following the second peak, the bolometric light curve 
declines at a rate of 20.011 mag day ! for about 300 days and 
then flattens for about 100 days before resuming a declining 
behavior. In addition to the overall fading we also find that the 
emission becomes redder with time, from initial colors soon 
after the second peak of g —rzr-—isz 0.1 mag, to 
g — r 7a 0.35 andr — i ~ 0.8 mag one year after the second 
peak. As we show below, this color evolution is primarily 
driven by an initially blue continuum, which eventually fades, 
leaving behind mainly strong nebular emission lines of calcium 
and oxygen. Our last data points have a magnitude 23.5 mag, 
deeper by about 0.5 mag than the pre-existing PS1/37 upper 
limits at the SN position. Since the SN is still detectable from 
the ground 2.5 yr after discovery, future ground and space 
observations will reveal the continued light-curve behavior, and 
potentially any underlying faint host. 

In Figure 2 we compare the light curves of SN 2016iet to 
those of both normal and unusually long-lived SNe I. 
SN 2016iet exhibits a distinct evolution from any of these 
events. Even its first peak alone is wider than the broad-lined 
Type Ic SN 1998bw (Galama et al. 1998) or even the longest- 
duration SN Ibn detected to date, OGLE-2014-SN-131 
(Karamehmetoglu et al. 2017). Other events match a similar 
overall timescale (e.g., the SLSNe SN 2015bn and PS1-14bj; 
Lunnan et al. 2016; Nicholl et al. 2016a, 2016b, 2018), but are 
either much less luminous or more rapidly fading from the 
initial brighter peak. The closest light-curve analogues are the 
Type Ic events iPTF15dtg (Taddia et al. 2019) and SN 2010mb 
(Ben-Ami et al. 2014), but even those do not share the double- 
peaked structure of SN2016iet. The SLSNe LSQl4bdq 
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Figure 2. Light curve of SN 2016iet in comparison to other types of hydrogen- 
poor and/or helium-poor SNe: SN 2010mb (Ben-Ami et al. 2014), iPTF15esb 
(Yan et al. 2017), iPTF1Sdtg (Taddia et al. 2019), SN 2015bn (Nicholl 
et al. 2016a, 2016b, 2018), PS1-14bj (Lunnan et al. 2016), SN 1998bw 
(Galama et al. 1998), and OGLE-2014-SN-131 (Karamehmetoglu et al. 2017). 
We also include a comparison to the He80 and Hel00 PISN models from 
Kasen et al. (2011), and to the hydrogen-rich PPISN candidate iPTF14hls 
(Arcavi et al. 2017). The comparison light curves have been smoothed for 
clarity, and shifted in phase to match the first peak of SN 2016iet. 


(Nicholl et al. 2015) and ¡PTF15esb (Yan et al. 2017) exhibit 
undulations in their light curves, but those are much less 
pronounced and occur on a much shorter timescale (20 days) 
than the double peaks in SN 2016iet. The Type IP SN 
iPTF14hls (Arcavi et al. 2017) has a light curve of longer 
duration than SN 2016iet and is the only known SN with 
undulations on timescales comparable to SN 2016iet; although 
these undulations could be caused by the same physical 
process, iPTF14hls is a hydrogen-rich event and therefore 
distinct from SN 2016iet. To conclude, we do not find a clear 
analog to SN2016iet in the existing SN sample with 
comparable behavior, brightness, and duration. 

In Figure 2 we also include a comparison to the PISN light- 
curve models for 80 and 100 Mo bare helium cores from Kasen 
et al. (2011). We find the peak brightness of SN 2016iet is 


Gomez et al. 


. mag x 1.15 (t/100d) 
; mag x 0.19 (t/100d) 


logio(Lpo;/erg s^!) 


6 x 10? 


Temperature / 1000K 


N 
o 


H 
o 
B 1.5 
o 
[m 
31.0 
© 
a 
m 


d 
Ui 


100 200 300 400 500 
Phase [rest days] 


Figure 3. Bolometric light curve (top), blackbody temperature (middle), and 
photospheric radius (bottom) of SN 2016iet. In all panels the dashed black line 
shows the result of a simultaneous fit to all epochs that contain multi-band data 
assuming a linear decline in temperature and radius. The green points are the 
result of interpolating the photometry in all bands to a common grid. We do not 
include the first peak due to the lack of color information. 


bounded by these two models, but the models lack the double- 
peaked structure and fade much more rapidly. 


3.2. Bolometric Light Curve 


We calculate the bolometric light curve of SN 2016iet, its 
photospheric radius, and temperature evolution by modeling the 
spectral energy distribution (SED) in each epoch by fitting a 
blackbody function using the Superbol code (Nicholl 2018). 
We further extrapolate the blackbody SED to account for the 
missing coverage in the UV and NIR. At phases earlier than 
75 days we assume a constant color to account for missing multi- 
band photometry. The resulting bolometric light curve, photo- 
spheric radius, and temperature evolution are shown in Figure 3. 

For the second light-curve peak we find a peak bolometric 
luminosity of =1.2 x 10” ergs !. The integrated radiated 
energy over the 650 rest-frame days following the second peak 
is E ® 2 X 10% erg. This does not include the first peak of 
the light curve for which we have no color information. 

We have gri-band photometry for phases between about 125 
and 420 days. At the time of the second peak we measure the 
photospheric temperature to be 7 = 9400 + 500K, and the 
radius to be R = (1.4 + 0.1) x 10!5 cm; we cannot provide a 
robust measurement of temperature before the second peak due 
to a lack of color information. Fitting all the data after the 
second peak simultaneously with an assumed linear decline, we 
find that the temperature decreases at a rate of 211 K day ' 
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Figure 4. Optical spectra of SN 2016iet spanning a rest-frame phase of 132-772 days relative to the first peak. The spectra are corrected for Galactic extinction and 
shifted to the rest frame of SN 2016iet using z = 0.0676. In some epochs we binned the spectra for clarity; the unbinned spectra are shown as shaded regions. The 
dashed horizontal lines indicate the zero level of each spectrum. The dominant lines in the spectra are marked. The data used to create this figure are available. 


and the radius decreases with a velocity of 2:300 kms ', 


shown in Figure 3. The bolometric light curve follows a decline 
of 0.0115 + 0.0002 mag day ! until day 420, and subse- 
quently flattens out to a decline rate of 0.0019 + 0.0018 mag 
day ! at a phase of 420-580 days. 

The data during the first peak are limited to only the Gaia 
G-band. To assess its bolometric luminosity and radiated energy 
we explore two extreme cases. First, assuming a constant 
temperature of about 8400 K during the first peak (the earliest 
temperature measurement at 84 days), we find a photospheric 
radius of =1.8 x 10P cm and a peak luminosity of =1.2 x 
10% ergs '. Second, assuming a radius equal to the smallest 
photospheric radius (7 x 10'*cm at 420 days) we find a 
temperature of 24,000 K and a peak bolometric luminosity of 
21.3 x 10% erg s !. This later luminosity corresponds to a total 
radiated energy during the first peak of Eaa ~ 10% erg, for an 
estimated rise time of 30 days. 


3.3. Spectral Features 


Our spectra of SN 2016iet span 132—772 rest-frame days and 
are shown in Figure 4. The spectra earlier than about 330 days 
are dominated by intermediate-width emission lines of Call, 
Or, [Carr], [O 1], [O rr], [O m], Mg], and a series of blended 
Fe lines near 5000 A. A strong blue continuum is also present. 
Beyond about 380 days, the blue continuum fades, although it 
does not completely disappear, and the spectra are dominated 
by strong emission lines of calcium, and oxygen, and weaker 
lines of sulfur and magnesium, leading to the color evolution 
seen in the light curves. 


The spectra show no clear evidence for hydrogen or helium 
associated with the SN, although we detect a feature at the 
location of Ha (11e confidence level) exclusively in the 
spectrum at 433 days. The origin of this feature is uncertain due 
to its faintness and presence in only a single spectrum. It could 
be due to ejecta interaction with a thin layer of hydrogen, or not 
be associated to Ho at all, but instead be produced by either 
[N II] A46548,6583, or Co II 46571, which has been previously 
observed in SNe Ia (Branch et al. 2004). 

In our final spectrum at a phase of 772 days we instead detect 
a narrow line at the location of Ha, shown in Figure 5. This line 
is unresolved (v < 230 km s^ P), and most likely associated with 
underlying star formation at the location of the SN. We measure 
its luminosity to be (3.9 + 1.5) x 1057 erg s !, which corre- 
sponds to a star formation rate of SFR ~ 3 x 10 ^ Ma yr ! 
(Kennicutt 1998). This value is close to the lower limit for dwarf 
galaxies, but similar to HII regions in the Small Magellanic 
Cloud (SMC; Kennicutt 1998; Hony et al. 2015). In the 
preceding spectrum (phase of 433 days) the SN was a few times 
brighter in the continuum, explaining the absence of this line in 
earlier spectra. 

The O I A7774 line is present in all the spectra and exhibits a 
broad red shoulder that is most pronounced in the early epochs, 
potentially due to contamination from Mg II AA7877, 7896 
(e.g., Nicholl et al. 2019). The emission lines in the complex 
feature around 4000 A are likely due to Call 143934, 3968. 
We see an unusual absorption doublet at -3850 A where the 
lines have a velocity width of 3600 + 550 km s": this feature 
could be caused by a P-Cygni profile of the Ca II doublet, but 
the exact nature of the line remains uncertain. 
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Figure 5. Zoomed-in on the location of Ha in the last two spectra of 
SN 2016iet. The spectrum at a phase of 433 days shows a broad line at the 
location of Ha. We fit this feature, and surrounding region, with a double 
Gaussian component plus a flat continuum and find a velocity width of 
2350 km s”*. At a phase of 772 days we detect narrow unresolved emission at 
the location of Ha likely due to star formation at the location of the supernova. 
We fit this line with a single Gaussian and find a velocity width <230 km s””. 


To measure the velocity width and integrated flux of the 
various lines we fit single or double Gaussian profiles to single 
and doublet lines, respectively; see Figure 6. We use the 
following lines: [Ca I] AA7291, 7324, Call 143934, 3968, OI 
A7774, [OI] 46300, [O n] A3727, and [O m] 144959, 5007, 
and find an average velocity of 3400 + 900kms !, with no 
evidence for significant change with time (Figure 7). The 
forbidden oxygen lines are all redshifted by =2800kms”' 
relative to the permitted oxygen lines and both the permitted 
and forbidden calcium lines (bottom panel of Figure 6). We 
present a possible interpretation for this shift in Section 3.5. 

Beginning from our first spectrum at a phase of 132 days, the 
lines of Ca II and OI increase in flux by a factor of about 2-3 
during the first 30 days, and then decline for the next 200 days 
by a factor of ~10. The [O Il] and [O m] line fluxes remain 
constant for the first 50 days, and then [O m] declines by a factor 
of ^15 for the next 250 days. The [OI] line flux appears to 
decrease monotonically by a factor of 10 for 200 days, starting 
with our first observation, and then remains constant for 100 
days. The strongest line is O I for most of our observations, until 
day 400, when it is overtaken by the [CaIAA7291, 7324 
doublet. The flux of this [Ca II] line remains constant within a 
factor of 2 throughout, until the last spectrum at phase 772 days. 

We measure the flux ratio of the [Ca 11] A47291,7324 to [O I] 
AA6300,6364 doublets at a phase of 433 days to be 
Licam/Lioy = 3.8 + 1.0. This is higher than for SNe IIb/ 
Ib/Ic, which have typical ratios of <2 (Milisavljevic et al. 
2017). This high ratio is also uncommon for SLSNe I, although 
comparable values have been measured for LSQ14an and 
PTF10hgi at late times (Nicholl et al. 2019). However, in those 
SLSNe the [Call] AA7291,7324 doublet was argued to be 
actually dominated by [Or] AA7319,7330 emission (Jerk- 
strand 2017); we do not believe this to be the case for 
SN 2016iet since we would expect this [O II] doublet to be 
redshifted like all of the other forbidden oxygen lines 
(Figure 6). High Licam/Lion ratios have been used as an 
indication of low progenitor masses, whereas we infer a high 
progenitor mass for SN 2016iet (Section 4). We caution that 
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Figure 6. Zoomed-in on the most prominent emission lines in the spectra of 
SN 2016iet. In some cases we exclude phases that lack spectral coverage or 
have a poor signal-to-noise ratio at the location of specific lines. Top: emission 
lines of Ca II, [Ca 1], and O I, which define the rest frame of the SN. Bottom: 
forbidden oxygen lines ([O I], [O IT], and [O M]), which are observed to be 
systematically redshifted by +2800 km s^! with respect to the SN rest frame. 


this relation has many caveats that depend on the phase of the 
observation and velocity of the ejecta, and holds only if [Ca I] 
and [OI] are produced in the same region, which is not 
necessarily the case for SN 2016iet (Fransson & Cheva- 
lier 1989; Jerkstrand 2017; Nicholl et al. 2019). 

The large ratio of Licam/Lro n is instead reminiscent of the 
calcium-rich gap transients (Kasliwal et al. 2012), which have 
values of =2.5-11 (Lunnan et al. 2017; Milisavljevic et al. 
2017). However, these transients are significantly less luminous, 
and evolve on much faster timescales, reaching peak magnitudes 
of only about —16 within 15 days after explosion, and entering 
the nebular phase in $100 days. Incidentally, calcium-rich gap 
transients have large projected offsets from their host galaxies 
(Lunnan et al. 2017), much like SN 2016iet (Section 5), but 
unlike SN 2016iet they are typically found in massive galaxies 
(Lunnan et al. 2017; De et al. 2018). 


3.4. Spectral Comparisons 


In Figure 8 we compare the spectra of SN 2016iet to those of 
various SNe I. The earliest spectra of SN 2016iet, at a phase of 
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Figure 7. Top: velocity width of the most prominent calcium and oxygen lines. 
No obvious evolution in the line width is observed over a period of nearly 2 yr. 
The solid and dashed horizontal lines mark the mean and standard deviation of 
the velocity width, 3400 + 900 km s_!. Bottom: fluxes of the same calcium 
and oxygen lines. The dashed color lines indicate the forbidden oxygen lines, 
which are redshifted by 2800 km s" from the rest of the lines. The velocities 
and fluxes are calculated from Gaussian fits to the line profiles shown in 
Figure 6. 


130-150 days exhibit emission lines of mainly calcium and 
oxygen, similar to the Type Ic SN 2007gr at a phase of 33 days 
(Valenti et al. 2008), although the emission lines in SN 2007gr 
are broader, with velocities of ~5000 kms '. SN 2007gr also 
does not exhibit the strong blue continuum present in 
SN 2016iet. SN 2016iet shows a much slower spectral evol- 
ution than typical SNe Ic, reaching the nebular phase only at 
about 400 days. The spectrum of SN 2016iet at a phase of 
350 days is similar to the 200 day spectrum of SN 2007Y, an SN 
Ib with high-velocity Ca II lines that was suggested to be powered 
by circumstellar medium (CSM) interaction (Stritzinger et al. 
2009). The spectrum of SN 2016iet at 462 days also resembles 
that of calcium-rich gap transients, for example SN 2005E (Perets 
et al. 2010; Lunnan et al. 2017) but at a phase of only about 60 
days, further emphasizing the fast evolution of calcium-rich gap 
transients as compared to SN 2016iet. We also show a comparison 
to PTF10hgi (Inserra et al. 2013; Quimby et al. 2018), an SLSN 
with one of the highest calcium-to-oxygen ratios (although as 
indicated above this ratio might be misinterpreted due to a 
possibly dominant [O II] line). 

The Type IIP SN iPTF14hls has a light curve with similar 
peak brightness and undulations to SN 2016iet, yet their spectra 
are completely different, where iPTF14hls is dominated by 
strong hydrogen lines (Arcavi et al. 2017). The Type Ic SN 
iPTFi5dtg also has a light curve similar in duration to 
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Figure 8. Comparison of three representative spectra of SN 2016iet (colored 
lines), including the early photospheric phase with a strong blue continuum 
(purple), the transition to the nebular phase at 328 days when the continuum 
begins to fade (pink), and the late nebular phase at 433 days with strong 
calcium and oxygen emission lines (yellow). We compare the evolution of 
SN 2016iet to other hydrogen- and/or helium-poor SNe (black). Strong telluric 
lines are marked with ©. The vertical lines mark the locations of dominant lines 
in the spectra of SN 2016iet. 


SN 201 6iet, but a spectrum unlike that of SN 2016iet that is 
instead similar to typical SNe Ic (Taddia et al. 2019). 

We show SN 2007bi as an example of a slowly evolving 
SLSN (Gal-Yam et al. 2009; Nicholl et al. 2013). The nebular 
spectra of SN 2007bi exhibit strong features of Fe-group 
elements, which we do not observe in SN2016iet. In 
SN 2016iet we observe much narrower lines than in SN 2007bi, 
closer to the velocities of ~2000-4000 km s! expected from a 
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PISN (Woosley 2017; Mazzali et al. 2019). We also plot a 
PISN model of a 100 Mẹ He core from Mazzali et al. (2019). 
However, the model differs significantly from the spectra of 
SN 2016iet, mainly due to a series of narrow iron lines. 

SN 2016iet is comparable only to SN 2010mb, with the 
caveat that SN 2010mb shows narrow emission lines of [O m] 
AA4959, 5007, [O u] A3727, and [OI] A5577 with a velocity 
with of 800 kms! (Ben-Ami et al. 2014). SN 2010mb also 
showed emission lines of calcium, oxygen, and magnesium 
with similar velocities to SN 2016iet, and a blue continuum at 
early phases, although not as prominent as in SN 2016iet (Ben- 
Ami et al. 2014). Unfortunately, the last spectrum of 
SN 2010mb was at only 131 days after explosion, preventing 
a direct comparison to the nebular spectra of SN 2016iet. 


3.5. Redshifted Oxygen Lines: A Light Echo? 


We find that all of the forbidden oxygen lines are redshifted 
by 2800 kms”! relative to the permitted oxygen lines and the 
permitted and forbidden calcium lines (Figure 6). This redshift is 
comparable to the velocity width of the lines. The redshifted 
lines are also symmetric around the redshifted velocity, rather 
than exhibiting a skewed velocity profile. The critical densities 
necessary to observe forbidden calcium and oxygen lines vary 
only by a factor of ~2, with n = 3.8 x 106 cm ? for [Or] 
26300,6364, and n = 8.8 x 10° cm ? for [Ca 1] A47291,7324 
(Jerkstrand 2017). This mild difference is unlikely to account for 
the different kinematics of the forbidden oxygen and calcium 
lines. 

Instead it appears that the permitted and forbidden oxygen 
lines must be arising from physically distinct regions, with 
the latter showing a lower density and an apparent recession. 
One possibility is that the forbidden lines are in fact an echo of 
the SN light from a circumstellar shell not yet reached by the 
ejecta, expanding at the observed redshift of 2800 km s !; such 
a velocity could arise in a giant eruption such as from the 
pulsational pair instability. If this shell was ejected before the 
SN, that would account for its lower density. Recently, Lunnan 
et al. (2018) detected an echo in the Mg 1112800 resonance line 
from a CSM shell around the SLSN iPTF16eh. In that case, the 
line showed a temporal evolution from blueshifted to red- 
shifted, due to the difference in light travel time to the observer 
between the front and back of the shell. In the case of 
SN2016iet, the lines are at all times redshifted, indicating that if 
they arise from an echo, we have missed the blueshifted phase. 
This sets constraints on the distance to the shell: to already see 
photons reflected from the far side of the shell in our earliest 
spectrum at 132 days, the inner radius of the shell must be 
located <1.7 x 10'’cm from the star. For a velocity of 
2800 km s_', this implies an ejection <20 yr before explosion. 

More difficult for this scenario is that the line remains visible 
over 300 days, which for a sharp injection of luminosity would 
require a very thick shell with a width (AR) several times 
greater than its inner radius. However, this may instead be a 
function of the fact that the light curve of SN 2016iet remains 
relatively flat for ~200 days. Given the long duration of the 
SN, it is then surprising that we do not see a blueshifted 
component. The unexplained absorption feature at 3850 A may 
offer a clue—if this is from blueshifted Ca II H and K lines, the 
implied velocity of the absorbing material is ~5000kms |, 
distinct from the redshifted forbidden oxygen lines. These lines 
could be coming from a distinct shell, or in fact not be related 
to the Call H & K lines. We defer further discussion of these 
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Table 1 

MOSFIT Parameter Definitions 
Parameter Definition 
Mo; Ejecta mass 
My; Radioactive nickel mass 
fxi Nickel mass as fraction of the ejecta mass 
Li Engine power for fallback model at 1 s 
ton Fallack power input constant up to this time 
Vej Ejecta velocity 
Ex Ejecta kinetic energy 
McsM CSM mass 
Ro CSM inner radius 
p CSM density 
Mys Neutron star mass 
Pspin Magnetar spin 
B, Magnetar magnetic field strength 
lexp Explosion time relative to first data point 
Tin Photosphere temperature floor 
MH host Column density in the host galaxy 
Ay host Extinction in the host galaxy 
K Optical opacity 
Ky Gamma-ray opacity 
c Uncertainty required for x =1 


lines until future detailed modeling, which is beyond the scope 
of this paper. 

In summary, SN 201 6iet shows a significant blue continuum, 
constant intermediate-width emission lines of calcium and 
oxygen with velocities of 3400 + 900 km sa high ratio of 
calcium to oxygen, no signs of hydrogen or helium, a peculiar 
set of redshifted forbidden oxygen lines that might be the result 
of a light echo, and a slowly evolving spectrum dissimilar to 
any normal SN Ic. 


4. Light-curve Modeling 


Due to the unusual photometric and spectrosopic nature of 
SN 2016iet, we explore a wide range of physical models to 
explain its origin, including radioactive decay, a magnetar and 
fallback accretion central engine, and circumstellar interaction. 
We fit the multi-band light curves using the Modular Open- 
Source Fitter for Transients (MOSFiT) Python package, a 
flexible Markov chain Monte Carlo (MCMC) code designed to 
model the light curves of transients using a variety of different 
power sources (Guillochon et al. 2018). We run each MCMC 
using an implementation of the emcee sampler (Foreman- 
Mackey et al. 2013). We test for convergence by ensuring that 
the models reach a potential scale reduction factor of «1.2 
(Gelman & Rubin 1992), which corresponds to about 
20,000-30,000 steps with 200 walkers, depending on the 
specific model. The parameters of the various models are 
summarized in Table 1. We note that MOSFiT implements 
simplified analytical one-zone models; the uncertainties 
presented in the following sections represent only the statistical 
errors on the fit. 


4.1. Radioactive Decay 


We begin with a radioactive nickel-cobalt decay model, 
which is the dominant power source in normal SNe I, although 
we note that this model cannot strictly account for the double- 
peaked structure of SN 2016iet. The input luminosity in this 
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Figure 9. Light curve of SN 2016iet with radioactive decay models for the full 
light curve (left) and only the first peak (right). Each line is a sample realization 
of the most likely models generated from MOSFiT. The resulting parameters 
for both models are provided in Table 2, and the models are discussed in 
Section 4.1. 


model is described by 
L, = Myi (enie "S + egge TTT, (1) 


where t is time, Mw; is the initial mass of 56Nj, eni = 6.45 x 
10” ergs™' Mz! and eco = 1.45 x 10% erg s! Ms! are heat- 
ing rates for ?9Ni and Co, respectively, and Tyni = 8.8 days 
and rc, = 111.3 days are the half-lives of “Ni and °Co, 
respectively (Nadyozhin 1994). 

The most likely models are shown in Figure 9. The models 
exhibit a broad light curve around the time of the first observed 
peak. But, as noted above, the model cannot account for the 
double-peaked structure, or for the flattening beyond 420 days. 
By the time of our latest data points the model underpredicts 
the observed brightness by about 1.5 mag. 

For this model we find M¿¡=64Mo, My; ® 3.7 Mo, 
Vej > 3500 km s^! (where the prior of the velocity was set by 
the observed width of the emission lines), and Ex = 4.4 x 
10% erg; the full list of parameters and their associated 
uncertainties is provided in Table 2. The model has a nickel 
fraction of fni = 0.06 + 0.02, consistent with broad-lined SNe 
Ic, which have typical values of fn; ~ 0.05-0.15 (Taddia et al. 
2019). However, this contrasts with the lack of evidence for a 
large fraction of Fe-group elements in the spectra of 
SN 2016iet, and the presence of a strong blue continuum. We 
therefore consider radioactive decay an unlikely model for the 
full light curve of SN 2016iet. 

Since we have no spectroscopic or color constraints during the 
first peak, and its duration is more comparable to normal SNe I 
(although still broader than even SN 1998bw, Figure 2), we 
explore a model in which only this portion of the data is explained 
by radioactive decay. In this case, we do not impose a limiting 
prior on Yai due to the lack of spectroscopic data. The best fit is 
shown in the right panel of Figure 9. We find Maj ~ 6 Mo, 
Msi © L3 Mc, (fyi © 0.24), vej > 2.0 x 10^ kms, znd Ex = 
14 x 10° ergs ! (Table 2). In this case both fw; and Ex are 
much larger than for SNe Ic (Taddia et al. 2019). We therefore 
consider this model unlikely, although it cannot be strictly ruled 
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Figure 10. Light curve of SN 2016iet with a magnetar (top) and fallback 
accretion (bottom) central engine model. Each line is a sample realization of the 
most likely models generated from MOSFiT. The resulting parameters are 
provided in Table 3, and the models are discussed in Sections 4.2 and 4.3. 


out. Naturally, this model still requires a distinct energy source for 
the remainder of the light curve. 


4.2. Magnetar Engine 


Here we explore the possibility that the long-duration light 
curve is powered by the spin-down energy of a millisecond 
magnetar (Kasen & Bildsten 2010; Woosley 2010). This model 
has been successful at explaining the diverse light curves of 
SLSNe I (see Nicholl et al. 2017 for a detailed description of 
the model and choice of priors). We note that this model cannot 
account for the double-peaked structure of the light curve, but 
modest modulations have been previously observed in SLSNe 
(Leloudas et al. 2012; Nicholl et al. 2015; Vreeswijk et al. 
2017; Blanchard et al. 2018). 

The most likely magnetar models are shown in the top panel of 
Figure 10, and the model parameters are listed in Table 3, where 
we use the observed width of the emission lines as a prior for 
the velocity of the ejecta. The model requires an initial rapid 
spin period of P ~ 1.5 ms, a magnetic field of B, ~ 3 x 10!* G, 
and a kinetic energy of Ex c 8 x 10° erg. These values are 
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Table 2 
Best-fit Parameters for the Radioactive Decay Models 
Parameter Prior Full Light Curve First Peak Units 
Mi; [1, 260] 64437 51:16 Mo 
Ni [0.001, 1] 0.06 + 0.02 0.2410? 
My; S 3:7203 1,5105 Ms 
Vo 3400 + 900° 3474* $55 19675+4984 kms 
Ej e 44738 13.7733 10% erg s 
tp [0, 500] 11848 17.9753 days 
Tas [1000, 15000] 76387383 37585071 K 
log(nu host) [16, 23] 17.0 + 0.07 17.9 + 1.7 cm? 
Av host? e] «0.002 «0.04 mag 
K [0.1, 0.2] 0.15 + 0.03 0.14 + 0.03 cmg! 
log(&,) [-4, 4] 1.6*1$ 1.5 + 1.9 cm?g ! 
logo [-4, 2] —0.53 + 0.03 —2.8 t 1.0 
Notes. Best model parameters, prior ranges, and lo error bars for the realizations shown in Figure 9. See Table 1 for parameter definitions. 
* These parameters were not fit for, but were calculated using all the posterior distribution samples of the fitted parameters. 
? Prior only used for the full light curve model; the prior for the first peak model was set to 1000-25000 km s~’. 
Table 3 
Best-fit Parameters for the Magnetar and Fallback Central Engine Models 
Parameter Prior Magnetar Fallback Units 
M [1, 260] 90% 93433 Mo 
log(L1) [51, 57] 552703 ergs! 
log(ton) [-4, 2] EE —2.1 € 1.3 days 
Mus [1.4, 2.2] 1.9+92 Mo 
Pis [0.7, 20] 1.5 + 0.4 ms 
Bi [0.01, 10] 3.2303 m 10^ G 
Ve; 3400 + 900 36991432 378378 kms”! 
Ef 7.5533 8.1482 10?! erg s 
feds [0, 500] 56.97 865 66.1715 days 
Tin [1000, 15000] 7648137 78213358 K 
10g(77.host) [16, 23] 18.6 + 1.8 18.5 + 1.7 cm ? 
Ay, host" <0.1 <0.08 mag 
K [0.1, 0.2] 0.16 + 0.03 0.14 + 0.04 e 
log(&,) [-4, 4] 1.5 + 1.7 —2.3 + 0.15 cm?)g ! 
logo [-4, 2] —0.75 + 0.04 —0.75 + 0.03 


Notes. Best model parameters, prior ranges, and lo error bars for the realizations shown in Figure 10. See Table 1 for parameter definitions. 
à These parameters were not fit for, but were calculated using all the posterior distribution samples of the fitted parameters. 


within the typical range inferred for SLSNe: P ~ 1-5 ms, B, ~ 
(0.1-3) x 10'%G, and Ex = (1-30) x 10° ergs’ (Nicholl et al. 
2017). However, this model requires a large ejecta mass of 
&:90 Mo, much higher than the typical values of 1-20 M. for 
SLSNe (Nicholl et al. 2017). We note that the magnetar model also 
nicely accounts for the late-time flattening of the light curves but, 
as mentioned above, does not strictly account for the double- 
peaked light-curve structure. 

A key challenge to the magnetar model, however, is that the 
inferred large ejecta mass falls in the regime for PISNe (i.e., 
helium core mass of 60-137 Mo; Heger et al. 2002), which is 
inconsistent with the required neutron star remnant. About 196 
of the posterior probability on Mej resides at masses below the 
PISN regime (i.e., in the PPISN regime) but in this case too the 
results are inconsistent since the star is expected to shed about 
half of its mass in the pulsational phase prior to explosion 
(Woosley 2017). 


10 


4.3. Fallback Accretion 


As an alternative central engine, we model the light curve with 
a fallback accretion model, in which some of the explosion debris 
remains bound and accretes onto the compact remnant (Dexter & 
Kasen 2013). The MOSFiT implementation of this model is 
described in Moriya et al. (2018). Our model is parameterized in a 
way that takes the engine power input L, (the luminosity at 1 s, 
assuming a power-law extrapolation), and calculates a constant 
input luminosity Lon = Li is up to time fọ» after which 
the input luminosity declines at a rate L c 1/3. The best-fit 
model for SN 2016iet has Lı = 10%? ergs and fon = 6905, 
which corresponds to an initial input luminosity of Lon =3 x 
10% erg sl during the first 690 s. We can estimate the total 
fallback mass from the accretion rate as 


ton t 
My = =| Í L.J + ‘rsa (2) 
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Assuming an efficiency of e = 0.1, and integrating for t= 
760 days, this yields Mg, = 1.7 Mo. 

The most likely fallback accretion models are shown in the 
bottom panel of Figure 10, and the model parameters are listed 
in Table 3. Overall, since the energy input rate is similar to that 
of a magnetar engine (i.e., t ^ /3 versus t 2), the resulting light 
curves and model parameters are nearly identical for the two 
models. In particular, we find Mej ~ 93 Mo. This model can 
also account for the late-time flattening of the light curves, but 
not for the double-peaked structure. 

This model therefore suffers from the same problems as the 
magnetar central engine model. Namely, the inferred progenitor 
mass lies in the PISN regime, expected to leave behind no 
remnant. Only about 6% of the posterior probability on Mej is 
in the PPISN regime, which still leads to an inconsistency, as 
discussed in the previous section. 


4.4. CSM Interaction 


We next fit the light curve with the CSM interaction model 
described in Chatzopoulos et al. (2013), where the radiation is 
produced by the interaction of the SN ejecta with a CSM that 
has either a wind (p cx r ?) or a shell (constant-density) profile. 
In this model, the double-peaked structure of the light curve 
could in principle be explained as being due to the forward and 
reverse shocks. We use the inferred radius of the photosphere 
(calculated in Section 3.2) to constrain the inner radius of the 
CSM, Ro, such that the radius of the photosphere is located 
inside the CSM at all times (i.e., Ro < Rpn). For this reason, we 
set the largest value of Ry to 54 au, the smallest robustly 
measured photospheric radius at a phase of 350 days. 

The results of the wind and shell interaction models are 
shown in the top and middle panels of Figure 11, and the full 
set of parameters is provided in Table 4. In both models we find 
a large CSM mass of Mcsm © 37 Mo. The ejecta masses are 
similarly large, My; ~ 85 and 45 Mo for the wind and shell 
models, respectively. The ejecta velocities of va ~ 4100 and 
2150 km s! in the wind and shell models are comparable to 
the observed line widths, although the value for the wind model 
is a better match (this is not influenced by the choice of prior). 
The inferred inner radius of the CSM is influenced by our 
choice of prior, with values of Ro 43 and «26 au, 
respectively. 

The CSM models are able to roughly capture the double- 
peaked structure of the light curve as the forward and reverse 
shocks. The models do not explain the late-time flattening, 
although this could be interpreted as the signature of an 
increase in the density of the CSM at larger radii relative to the 
simple density profiles assumed here. Finally, we note that the 
CSM interaction model is expected to produce a blue 
continuum, which we indeed observe in the spectra 
(Figure 4). Thus, the CSM models provide an adequate 
explanation for the light curves and spectra, and require a 
progenitor mass shortly before explosion (i.e., a CO core mass 
of Mej + Mcsm) of =85-120 Mo, squarely in the PISN regime. 


4.5. Shock Breakout plus CSM Interaction 


Finally, we explore a variant of the CSM interaction model 
in which the first peak is due to shock breakout from an inner 
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Figure 11. Light curve of SN 2016iet with CSM interaction models for a wind 
density profile (p ac r ?; top), a constant-density shell profile (middle), and a 
shell profile fit to only the second peak and beyond (bottom). Each line is a 
sample realization of the most likely models generated from MOSFiT. The 
resulting parameters are provided in Table 4, and the models are discussed in 
Sections 4.4 and 4.5. 
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Table 4 

Best-fit Parameters for the CSM Interaction Models 
Parameter Prior Wind CSM Shell CSM Shell CSM Excluding First Peak Units 
Mi; [1, 260] 84.6*183 45.41159 18.9+190 Mo 
Mesm [1, 260] 36.31% 37.9*125 33.0764 Mo 
Vi [1000, 10000] 41001595 21477398 2599+539 kms” 
Ef 873 12:04 0203 10?! erg s_! 
Ro [1, 54.2] 42.8182, 26.37 138 44.4 + 7.0 au 
log(p) [-14, -9] —11.1 + 0.2 —12.0 + 0.2 —11.8 4 0.2 gem 
ias [0, 500]? 60.071232 128*13 107.144 days 
Ti [1000, 15000] 77533416 78894373 7652+33 K 
log(11H, host) [16, 23] 18.2412 18.7474 18.6 + 1.8 em? 
AV, host’ p «0.06 «0.30 «0.13 mag 
k [0.1, 0.2] 0.14+00 0.14+90 0.13*005 e 
logo [-4, 2] —0.63 + 0.04 —0.61 + 0.04 —0.75 + 0.04 


Notes. Best model parameters, prior ranges, and 1c error bars for the realizations shown in Figure 11. See Table 1 for parameter definitions. The maximum on the 
prior of Ro is motivated by the fact that the photosphere radius should be greater than the inner radius of the CSM. 

* These parameters were not fit for, but were calculated using all the posterior distribution samples of the fitted parameters. 

P Priors only used for wind and shell CSM models. The priors for the “shell CSM excluding first peak” model were set to [34.65, 54.2] au and [113.7, 100.9] days, 


respectively. 


dense CSM and the second peak and subsequent light curve are 
powered by CSM interaction as in the previous section. 

We model the first peak with the shock breakout model of 
Chevalier & Irwin (2011, Equations (3) and (5) in that paper) from 
which we can infer the breakout radius Rj, as well as Mej, Mcsm, 
and Ex. The input parameters are the light-curve rise time fise, the 
radiated energy Eraa, the shock velocity vs, and an optical opacity 
x. We estimate tise ~ 30 days from a second-order polynomial fit 
to the first peak, which is similar to the 28 days separation 
between the first and brightest Gaia data points. We calculate 
Eaa ~ Lp X trise ~ 10% erg, where Lp ~ 1043-1044 erg s! 
(calculated in Section 3.2). We set the value of « to 0.1-0.2, 
typical values assumed for hydrogen-poor material. The resulting 
ejecta and CSM parameters as a function of vs, are shown in 
Figure 12. We find that self-consistent solutions require 
va, 2200-2800 km s~, comparable to the observed emission 
line widths in the spectra. The inferred breakout radius is 
Kua © 40-50 au, the ejecta mass is My ~ 15-28 Mo, and the 
CSM mass of the breakout region is =1-4M¿. We find a 
parameter D, = 8-16 (Chevalier & Irwin 2011), depending on 
the assumed opacity, corresponding to a mean density in the 
breakout region of p ~ 2 x 10 !! gem |. 

We model the rest of the light curve with interaction with a 
shell-like CSM as in Section 4.4, but with an inner radius 
Ro Z 35 au, the smallest inferred value of Ry, shown in 
Figure 13. We find Mcsm = 33 Mo, M, = 19 Me, and 
Vj ® 2600kms '; see the bottom panel of Figure 11 and 
Table 4 for the model light-curve and parameter values. 
Comparing these results to the breakout analysis of the first 
peak we find consistent values of vej, Mej, and Ex. Considering 
the overall simplified breakout model and the rough estimate of 
the rise time and radiated energy of the first peak, we consider 
the overall agreement to indicate that this is a plausible 
interpretation of the double-peaked structure. 

This result is also supported by the density inferred from the 
two components. In Figure 13 we plot the density profile from the 
breakout model, which, following Chevalier & Irwin (2011), we 
assume to be a wind profile, and the constant density inferred 
from the CSM interaction model of the second peak. We find that 
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the two join smoothly at Ry, Y Resmin ~ 40 au with a density of 
10712 g cm ^. 


4.6. Summary of the Light-curve Models 


We modeled the light curves of SN 2016iet with a broad 
range of energy sources. We find that a radioactive decay 
model fit to the entire light curve requires a large ejecta mass of 
64 Mo and a nickel mass of 3.7 Mo. This model has several 
disadvantages; it cannot account for the double-peaked 
structure or the late-time flattening of the light curves, it 
severely underestimates the latest measurements, and it requires 
a nickel mass fraction that is not supported by the spectroscopic 
data. We therefore consider this model unlikely. A radioactive 
model fit for just the first peak requires an unusually large 
nickel fraction (70.24), making it an unlikely explanation. 
While the latter possibility cannot be ruled out, it still requires 
an explanation for the rest of the light curve. 

The magnetar and fallback accretion central engine models 
both require a large ejecta mass of 2:90 Mo. These models do 
not explain the double-peaked light curve structure, but they 
nicely capture the late-time flattening. The inferred engine 
parameters for the magnetar central engine are similar to those 
of SLSNe, but the ejecta mass is much larger. A key theoretical 
problem for the engine models, however, is that progenitors 
with a helium core mass this high are expected to produce 
PISNe, which do not leave behind a neutron star or black hole 
remnant. 

In the ejecta-CSM interaction models, the double-peaked 
structure is accommodated by emission from the forward and 
reverse shocks. This model can also account for the blue 
continuum seen in the spectra in the first year. The resulting 
CSM mass is 2:37 Mo and the ejecta mass is Mej ~ 45-85 Mo, 
depending on the density profile of the CSM. Therefore, the 
progenitor had an inferred mass of 285-120 Mo shortly before 
explosion, after losing its hydrogen and helium layers. This 
corresponds to a progenitor ZAMS mass of ~150-260 Mo, the 
regime expected to produce PISNe. The late-time flattening is 
not captured in these models, but it can in principle be 
attributed to a change in the density profile at larger radii. 
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Figure 12. Solutions to the shock breakout equations of Chevalier & Irwin 
(2011) used to calculate the kinetic energy Ex,51, ejecta mass Mej, CSM mass 
Mcsm, and breakout radius Ry, from the inferred radiated energy E,aa, rise time 
fise, Opacity &, and shock velocity vs. The cyan band represents the solution 
with a fixed radiated energy of 10°° erg, and the purple band represents the 
solution with a fixed rise time of 30 days. The regions where both bands 
overlap are the allowed solutions. The green and blue boundaries are 
determined by the limits in the allowed opacity « = [0.1, 0.2]. The red data 
points are the result from the MOSFiT CSM interaction models excluding the 
first peak shown in the bottom panel of Figure 11. We find that the two fits are 
fully self-consistent. 


In the breakout+interaction model, we find an ejecta mass of 
=19 Me, an inner breakout region with a mass of =3 Mọ 
extending to about 40 au, and an outer CSM with a mass of 
=33 Mo. This model therefore requires the lowest progenitor 
mass, with a CO core of =55 Mo. This variation of the CSM 
model provides a consistent result in terms of the ejecta and 
shock velocity, ejecta mass, and kinetic energy from the 
breakout-powered first peak and the interaction-powered 
second peak (Figure 12). 

The implication of the CSM models (regardless of variant) is 
a mass of Mcsm ® 35 Mo located between inner and outer 
radii of 40-150 au (e.g., Figure 13). For an ejection velocity 
of =100kms”*, this requires a mean mass-loss rate of M ~ 
TMs vr" spanning the final =2—7 yr before explosion. If the 
ejection velocity was instead 1000 km s^! the corresponding 
timescale over which the CSM was ejected would be only 
70-260 days before explosion, with a mean rate of 
0.2 Mo day. 

The key result from our modeling is that, regardless of the 
model, the inferred total mass (ejecta or ejecta-CSM) is well 
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Figure 13. Inferred density profile of the CSM from the shock breakout model 
of Chevalier & Irwin (2011) for the first peak of the light curve (gray band), 
and a shell density profile CSM model for the second part of the light curve 
(red band; bottom panel of Figure 11). The two density profiles are fully 
consistent. 


outside the range inferred for normal SNe I, including the most 
energetic Type Ic events and SLSNe, which have Mej $ 20 Mo 
(Modjaz et al. 2016; Nicholl et al. 2017). We show the model 
posterior distributions of the total mass in Figure 14 in 
comparison to theoretical predictions for various explosion 
mechanisms. Stars with a helium core mass of 2-40 MG are 
expected to produce an ordinary core-collapse SN, or in some 
cases a failed SN. Stars with a helium core mass of 40-60 Mo 
are expected to produce a PPISN, while those with a helium 
core mass of 60-137 Mẹ are expected to produce a PISN 
(Woosley et al. 2007). All of the models for SN 2016iet place 
the progenitor in the PPISN or PISN regimes, especially the 
CSM models, which provide the most consistent explanations 
for the light curves and spectra. 


5. Host Galaxy 


No galaxy is detected at the location of SN 2016iet in archival 
PS1/37 images to 30 limits of g = 23.1 mag, r = 23.0 mag, 
i 222.7 mag, z 3 22.0 mag, and y = 21.0 mag. The nearest 
galaxy with a low probability of chance coincidence of Pec ® 
0.1, is located 11.1 + 0"3 from the SN position (Figure 15). 
We measure a redshift for this galaxy of z — 0.0671 + 0.0001 
from several prominent narrow emission lines (Ho, HØ, [O 1] 
AA3726,3729, and [O m] A44959,5007; see Figure 15). This 
represents a velocity difference of 140 + 85 km s! relative to 
the redshift inferred from the emission lines in SN 2016iet. We 
therefore conclude that SN 2016iet is associated with this galaxy. 

We calculate a projected physical separation of SN 2016iet 
to this host galaxy of 16.5 + 0.6 kpc. Relative to the half-light 
radius of the host galaxy, Re = 3.8 + 0.4 kpc, this corresponds 
to a normalized offset of D/R, — 4.3 + 0.4. Kelly & Kirshner 
(2012) found that SNe Ic tend to occur close to their host 
galaxies, with D/R, < 2.4. Similarly, SLSNe exhibit offsets of 
D/R, S 1.6 (Lunnan et al. 2015), while long gamma-ray bursts 
(LGRBs) have normalized offsets of D/R, X 1.8, where only 
one out of a sample of 105 has a separation D/R, ~ 4 
(Blanchard et al. 2016). Calcium-rich gap transients show 
comparable host offsets, with a median of D/R, = 10 (Lunnan 
et al. 2017), but these transients tend to occur in massive 
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Figure 14. Inferred ejecta or ejecta+CSM masses from the different model fits 

to the light curve. The shaded regions mark the theoretical ranges for core- 

collapse SNe (green), PPISNe (blue), and PISNe (pink). The posterior mass 


distributions for the models of SN 2016iet all overlap the PPISN or PISN 
ranges. 


galaxies (Lunnan et al. 2017; De et al. 2018), unlike the host of 
SN 2016iet. In general, massive progenitors tend to occur at 
small offsets from their hosts, whereas SN 2016iet clearly 
violates this trend. 

We extract the host's magnitudes by performing photometry 
on stacked PS1/37 grizy template images using a constant 
circular aperture, to ensure uniform sampling in all filters. We 
include additional ugr measurements from our IMACS and 
LDSS3c images. The resulting magnitudes are listed in 
Table 5, along with the corresponding absolute magnitudes. 
We find M, ~ —17.5 mag, or L, ~ 0.02 L. (Montero-Dorta & 
Prada 2009), identical to that of the SMC. If SN 201 6iet instead 
originated in an underlying satellite galaxy, then the brightness 
of this galaxy would be M, = —13.9 mag, or L, < 0.003 Lx. 

The flux ratio of Lya /Lyg = 2.3 + 0.6 is consistent with the 
expectation for case B recombination (Lga/Lgg = 2.86), 
indicating no evidence for significant host galaxy extinction. 
We infer the metallicity using the Kaa method (Kobulnicky 
et al. 1999), which is double-valued, but the lack of detectable 
[N II] emission and the low luminosity of the host galaxy point 
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Figure 15. Image of SN 2016iet and its most likely host galaxy taken with the 
LDSS3c instrument in i-band on 2018 July 9. The blue lines mark the position 
of the slit during two spectral observations that included the host. The resulting 
host spectra are shown in the lower panel. 


toward the lower branch solution. This gives a value of 
12 + log(O/H) = 7.807013, or Z= 0.144003 Z. (Asplund 
et al. 2009). Given the large offset, the metallicity at the 
location of SN 2016iet is likely even lower. We stress that this 
metallicity is low even in comparison to the host galaxies of 
SLSNe, which occur preferentially in low-metallicity galaxies, 
but mostly span 12 + log(O/H) = 7.8-8.6 (Lunnan et al. 
2014); see Figure 16. Similarly, only 596 of the SN Ic sample 
from Modjaz et al. (2019) have metallicities lower than the 
value inferred for SN2016iet, namely SNe 2010hie and 
2011rka. 

To extract additional information about the host we model its 
SED using Prospector, a software package designed to fit 
galaxy SEDs by taking into account dust attenuation and re- 
radiation, nebular emission, stellar metallicity, and a piecewise 
linear non-parametric star formation history (Leja et al. 2017). 
We model the SED with seven bins of star formation history, 
spaced evenly in logarithmic time from t = 0 to the age of the 
universe at z — 0.0671, allowing a two-component dust screen 
model with a flexible dust attenuation curve (Charlot & 
Fall 2000). The resulting SED and Prospector model are 
shown in Figure 17. We find a rising star formation history, 
where the most recent star formation rate is SER 39 uy; = 
0.497029 M... yr !. We find a mass of log(M/M.) = 8.56017, 
and a stellar metallicity of Z = 0.08*006 Z... consistent with 
the gas-phase metallicity we inferred from emission line fluxes. 
We also fit a parameter permitting ionizing radiation to escape 
from the galaxy, instead of powering nebular emission. 
Inclusion of this parameter notably improves the quality of 
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Table 5 
Host Galaxy Properties 
Value Units Value Units 

B; 0.12 M. —17.50 + 0.12 mag 
z 0.0671 + 0.0001 M, —17.35 + 0.34 mag 
D 11.1 + 0.3 arcsec Ey 0.28 + 0.04 Ly 
D 16.5 + 0.6 kpc E; 0.076 + 0.004 Ly 
50 2.57 + 0.26 arcsec L, 0.0234 + 0.0004 Ly 
D 140 + 85 km s_! L; 0.020 + 0.002 Ly 
u 20.27 + 0.15 mag (IMACS) L. 0.013 + 0.001 Ly 
2 19.96 + 0.06 mag (3PI) Lia 15.5 + 1.4 10% erg s^! 

19.94 + 0.14 mag (LDSS) Lug 6.8 + 1.8 10% erg s^! 
r 20.00 + 0.05 mag (3PI) Lio 11137263729 13.5 + 13 10% erg s~! 
r 20.04 + 0.02 mag (LDSS) Lio 1950 3.7 + 1.0 10% erg s~! 
i 19.96 + 0.09 mag (3PI) Lio 11115007 12.2 + 2.6 10% erg s^! 
i 20.11 + 0.14 mag (LDSS) 12 + log( O/H) 7.807013 
z 19.98 + 0.12 mag (3PI) Z3 0.147003 Zo 
y 20.12 + 0.34 mag (3PD SFR; 0.497019 Mo yr! 
M, —17.21 + 0.15 mag SFR; 0.087933 Mo yr! 
My —17.53 + 0.06 mag log(M) 8.56+91 Mo 
M, —1745 + 0.02 mag VA 0.08006 Zo 
M; —17.53 + 0.09 mag et 


Note. List of measured and derived properties of the host galaxy of SN 2016iet. P.o is the probability of chance coincidence, z is the redshift, D is the offset of 
SN 2016iet from the host center, rso is the half-light Petrosian radius in the r-band, v, is the relative radial velocity difference between SN 2016iet and the host, ugrizy 
are aperture magnitudes corrected for Galactic extinction, M,,¢,;, are the host absolute magnitudes, L. are the host luminosities relative to L in each band (Montero- 
Dorta & Prada 2009), L, are the integrated luminosities of several key emission lines (Ha, H, [O r1], [O m]), 12 + log(O/H), and Z>3 are values for the host 
metallicity inferred from the R23 formalism. SFR1 and SFR2 are the earliest 30 and 100 Myr of star formation from prospector models. The mass and metallicity 


of the galaxy log(M) and Z are also inferred from the prospector models. 
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Figure 16. Inferred metallicity of the host galaxy of SN 2016iet (green) relative 

to the cumulative distribution of the host galaxy metallicities for Type Ic SNe 

(Modjaz et al. 2019) and SLSNe (Lunnan et al. 2014). Also shown in red is the 

host metallicity for SN 2010mb (Ben-Ami et al. 2014). 


the fit to the photometry. We infer a high escape fraction 20.9 
from this method; we caution that this may be driven in part by 
deficiencies of the nebular emission model at low metallicities 
and high ionization states rather than true ionizing continuum 
leakage. 

In our final spectrum at a phase of 772 days we detect 
unresolved Ha emission at the location of SN 2016iet 
(Figure 5), most likely originating from star formation. The 
Ha luminosity corresponds to SFR ~ 3 x 10 ^ Mo yr’. This 
is much smaller than for dwarf galaxies (1.5 Mẹ yr ' for the 
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Figure 17. Distribution of the most likely Prospector models to the SED of 
the host galaxy of SN 2016iet. The blue data points are the observed 
photometry listed in Table 5. The shaded gray region shows the +10 scatter in 
the posterior of Prospector models. The inset shows the recovered star 
formation history. 


SMC; Rezaeikh et al. 2014), but similiar to H II regions in the 
SMC. For example, the NGC 346/N66 nebula, which has an 
SFR = 4 x 10? M. yr ! (Hony et al. 2015). Suggesting the 
progenitor of SN 2016iet could have formed in situ in an 
isolated H II region. Crowther (2013) finds that seven out of 10 
hydrogen-poor SNe in their sample are associated with an H II 
region. If a compact, isolated HII region is present at the 
location of the SN, a very high-mass progenitor (>75 M.) is 
favored. Since for high-mass isolated clusters, gas has already 
been removed, we would not expect an SN coincident with the 
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HII region, unless this was from a very massive progenitor 
(i.e., a lifetime shorter than the duty cycle of giant H II regions.) 

To summarize, SN 2016iet is associated with a low- 
metallicity dwarf galaxy, with an unusually large offset relative 
to other H/He-poor massive star explosions (SNe Ic, SLSNe, 
and LGRBs). Since the range of possible progenitor core 
masses for SN 2016iet is 55-120 Mo, which corresponds to 
ZAMS masses of “120-260 Mo, we conclude that the 
progenitor of SN 2016iet must have formed in situ. Stars of 
these masses have lifetimes of 3-4 Myr (estimated from 
stellar evolutionary tracks; Choi et al. 2016). If the large offset 
was due to ejection, this would require a velocity of 
24000 km s !, much higher than even the fastest hypervelocity 
stars in the Milky Way (1200 kms ^!; Boubert et al. 2018). 
This therefore indicates that the progenitor of SN 2016iet likely 
formed in situ, in a dim satellite galaxy or an isolated star 
cluster. This is further evidenced by our detection of narrow 
Ha at the location of the SN. 


6. SN 2016iet as a PPISN or PISN 


We infer that the progenitor of SN 2016iet was a CO core 
with a mass near the end of its life of 255 Mo and potentially 
up to 120 Me. This is the range of core masses that are 
expected to produce PISNe or PPISNe (Woosley et al. 2007). 
We further find that the environment of SN 2016iet has a 
metallicity of Z ~ 0.1 Zo, in agreement with models for PISN 
and PPISN progenitors that suggest a threshold of <0.1 Zo 
(Yusof et al. 2013; Owocki 2015; Yoshida et al. 2016). The 
large ratio of Licam/Lion = 3.8 + 1.0 we observe in 
SN 2016iet is also in agreement with models of PISNe, which 
show a ratio of ~2.5 for a 100 Mẹ He core model (Mazzali 
et al. 2019). Within the context of our models, CSM interaction 
provides the most self-consistent explanation for the data; the 
radioactive model is not consistent with the spectra, while the 
engine models are not consistent with the notion that PISNe do 
not leave remnants. 

Our minimal mass interaction model has an inferred CSM 
mass of about 35 Mo, of which about 33 MG are located in a 
shell that spans about 40-150 au from the progenitor, and a 
few Mo located at smaller radii (in the breakout region). The 
ejecta mass is 220 Mo, and the overall CO core mass is in the 
PPISN regime (Figure 14). The other interaction model variants 
have identical CSM masses and locations to the this model, but 
much larger ejecta masses of about 45-85 Mo that place the 
progenitor in the PISN regime (Figure 14). 

Light-curve models for PISNe from He cores, in which the 
emission is powered by radioactive decay and the diffusion of 
thermal energy (e.g., Kasen et al. 2011), do not match the data 
for SN2016iet. Kasen et al. (2011) suggest that CSM 
interaction may provide an additional energy source, but they 
do not model this effect or make specific predictions for the 
structure of the CSM (e.g., mass, density, ejection timescale). 
To match the inferred properties of SN 2016iet, stripped PISN 
progenitors would have to experience severe mass loss prior to 
the eventual explosion. 

Theoretical models of PPISNe have a wide range of 
predicted observational properties. Woosley (2017) carried 
out a detailed study of PPISNe from helium cores of 30-64 Mo 
and their pulsation histories. We compare SN 2016iet to the 
PPISN model of a 62 Mo helium core. In this model 28 MG of 
material are ejected in a pulsational mass-loss event about 
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2700 yr prior to core collapse, leaving behind a 34 Mo remnant 
that may experience the mass loss of an additional few Mẹ in 
the final years before explosion, and then most likely collapse 
to a black hole. While the large mass ejection and remaining 
core mass are comparable to what we infer for SN 2016iet, the 
mass ejection timescale is off by nearly three orders of 
magnitude. Moreover, our inference of a large ejecta mass for 
SN 2016iet would require the core-collapse explosion to be 
successful. 

Thus, while at face value SN 2016iet matches the predictions 
for PPISN/PISN explosions in terms of the inferred progenitor 
mass and metallicity, in detail, it appears to challenge these 
models by requiring that about half of the progenitor's mass be 
lost in the final years before explosion. This suggests that either 
SN 2016iet was formed by a completely different mechanism, 
or that PPISNe/PISNe progenitors may undergo different 
behavior than suggested by current theoretical models. 


7. Conclusions 


We present a detailed optical follow-up of SN 2016iet, 
including the determination of its redshift and identification of 
its host galaxy. We classify this event as a unique Type I SN 
due to the lack of H and He. The combination of its unusual 
light curve, spectra, and environment is unprecedented; the 
only potential analog in the existing literature is SN 2010mb, 
but even it does not share all of the same properties (e.g., 
luminosity, duration, and low-metallicity host). The key 
observational properties of SN 2016iet are as follows. 


1. The light curve exhibits a double-peaked structure with 
nearly equal brightness peaks of about —19 mag 
separated by about 100 rest-frame days. 

2. The subsequent decline in brightness is slow, =3 mag in 
=650 rest-frame days, with the SN remaining detectable 
at present. 

3. The spectra are dominated by intermediate-width 
(=3400 kms *) emission lines of calcium and oxygen, 
superposed on a strong blue continuum in the first year, 
and lack hydrogen or helium. 

4. The ratio of [Ca N] to OI of about 4 in the late (nebular) 
spectra is much larger than for normal stripped core- 
collapse SNe and SLSNe. It is reminiscent of Ca-rich 
transients, although the duration and luminosity of 
Ca-rich transients are far lower than those of SN 2016iet. 

5. The forbidden oxygen lines are redshifted by 
=2800 km s" with respect to the permitted oxygen lines 
and permitted and forbidden calcium lines. The observed 
redshift may be indicative of a light echo from a 
circumstellar shell not yet reached by the ejecta. 

6. The SN is located at a large offset (16.5 kpc, 4.3 effective 
radii) from a low-metallicity dwarf galaxy. This offset is 
substantially larger than for other types of stripped 
envelope explosions (Type Ib/c SNe, SLSNe, LGRBs). 

7. The host galaxy has a low metallicity of Z ~ 0.1 Zo, 
unusually low even in comparison to SLSN and LGRB 
host galaxies. 

8. The detection of a weak and narrow Ha line at late time 
(phase of 772 days) coincident with the SN location 
indicates underlying star formation activity, either from a 
dwarf galaxy with L < 0.003 L,, or an isolated star- 
forming region. 
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We explore several models for the energy source powering 
SN 2016iet—radioactive decay, central engine, and CSM 
interaction—and find that, regardless of the model, the inferred 
mass of the progenitor (CO core) prior to explosion is in the 
range of =55-120M;. This suggests a ZAMS mass of 
120-260 Mo, which along with the low metallicity, places 
the progenitor in the regime of PPISN/PISN explosions. This 
is the first PPISNe/PISN candidate with an inferred progenitor 
mass and metallicity in the range predicted by theoretical 
models. 

In detail, the most consistent model for the light curve and 
spectra is CSM interaction, in which the CSM has a mass of 
=35 Mə and is located over a range of about 40-150 au of the 
progenitor. This indicates that the CSM was ejected within a 
decade before the explosion at a rate of 2:7 Mẹ yr ?. In the 
model of breakout+interaction, an additional few Mẹ of 
material are required in even closer vicinity to the progenitor. 

While the mass of the progenitor and its inferred metallicity 
are overall consistent with the predictions of PPISN/PISN 
models, we find that, in detail, the inferred properties are 
challenging to explain in the framework of existing predictions. 
In particular, PPISN models in the range of masses inferred 
here can indeed eject about half of their CO core mass, but are 
predicted to do so thousands of years before explosion 
(Woosley 2017), not in less than a decade, as required for 
SN 2016iet. Additional mass loss may occur in the final years 
before the final core collapse, but at a level of only a few Mo 
(Woosley 2017; Aguilera-Dena et al. 2018), an order of 
magnitude too low for SN 2016iet. Finally, Woosley (2017) 
also argues that the remaining cores of 35-45 M¿ have a 
large binding energy that may make them difficult to explode, 
and lead instead to massive black holes; this is again in tension 
with the inferred ejecta mass of ~20-85 Mo for SN 2016iet. In 
the context of PISN models, CSM interaction has been 
mentioned as a potential energy source for explosions of He 
cores, but this has not been explored in detail (e.g., Kasen et al. 
2011). 

We end with the following cautionary note. SN 2016iet may 
indicate that existing predictions for PPISN/PISN explosions 
are flawed. Namely, either PPISN/PISN progenitors can eject 
about half of their mass in the final years before explosion 
(rather than several thousand years before explosion), or 
PPISNe can actually leave behind compact remnants in the 
form of a millisecond magnetar or an accreting compact object 
(as required in the context of the engine models for 
SN 2016iet). Alternatively, SN2016iet could instead be 
powered by a completely different mechanism not explored 
in this work, such as a binary merger involving at least one 
CO core. 
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Appendix 
Data Tables 


Here we present the optical photometry of SN2016iet in 
Table 6 and optical spectroscopy in Table 7. 
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Table 6 

Optical Photometry 
UT Date MJD g (err) r (err) i (err) 
(year month day) (mag) (mag) (mag) 
2016 May 19 575214 eU U 19.7 
2016 Jul 11 57580.8 Ut UD et 
2016 Nov 14 5706.6 
2016 Nov 14 57706.8 
2016 Dec 10 57732.0 
2016 Dec 10 57732.1 
2017 Jan 05 57758.3 
2017 Feb 1 57785.5 
2017 Feb 1 57785.5 
2017 Feb 1 57785.5 
2017 Feb 1 57785.5 
2017 Feb 23 57807.5 
2017 Feb 23 57807.5 
2017 Feb 23 57807.5 
2017 Feb 23 57807.5 et ve et 
2017 Mar 6 57818.4 et UD 18.75(0.02) 
2017 Mar 6 57818.4 sE Ut 18.75(0.02) 
2017 Mar 6 57818.5 et U 18.79(0.03) 
2017 Mar 6 57818.5 et Ut 18.77(0.02) 
2017 Mar 6 57818.5 ... Ut 18.77(0.03) 
2017 Mar 16 57828.5 et Ut et 
2017 Mar 16 57828.5 
2017 Mar 16 57828.5 
2017 Mar 16 57828.5 et Ut et 
2017 Mar 16 57828.5 sE ... 18.78(0.05) 
2017 Mar 16 57828.5 ... a 18.76(0.04) 
2017 Mar 16 57828.6 ... Ut 18.80(0.06) 
2017 Mar 16 57828.6 ... Ut 18.80(0.06) 
2017 Mar 25 57837.5 et UD Ut 
2017 Mar 25 57837.5 
2017 Mar 25 57837.5 
2017 Mar 25 57837.5 
2017 Mar 31 57843.5 
2017 Mar 31 57843.5 
2017 Mar 31 57843.5 
2017 Mar 31 57843.5 e UD Ut 
2017 Apr 15 57858.5 et Ut 18.91(0.06) 
2017 Apr 15 57858.5 Ut UD 18.84(0.06) 
2017 Apr 15 57858.5 et Ut 18.86(0.07) 
2017 Apr 15 57858.5 Ut UD 18.86(0.06) 
2017 Apr 15 57858.5 Ut Ut Ut 
2017 Apr 15 57858.5 
2017 Apr 15 57858.5 
2017 Apr 15 57858.5 
2017 Apr 21 57864.5 
2017 Apr 21 57864.5 
2017 Apr 21 57864.5 
2017 Apr 21 57864.5 
2017 May 21 57894.5 
2017 May 21 57894.5 
2017 May 21 57894.5 
2017 May 21 57894.5 eU UD Ut 
2017 May 23 57896.3 19.14(0.10) 19.07(0.09) 19.23(0.13) 
2017 May 24 57897.3 19.22(0.10) 19.10(0.10) 19.10(0.13) 
2017 May 28 57901.3 19.23(0.09) 19.17(0.09) 19.21(0.15) 
2017 May 29 57902.3 19.35(0.10) 19.28(0.13) 19.40(0.15) 
2017 Jun 06 57910.2 19.29(0.11) 19.24(0.14) 19.29(0.12) 
2017 Jun 06 57910.3 Ut Ut 19.37(0.15) 
2017 Jun 06 57910.3 et UD 19.38(0.15) 
2017 Jun 06 57910.3 eU Ut 19.28(0.15) 
2017 Jun 11 57915.2 19.31(0.06) 19.26(0.09) 19.27(0.13) 
2017 Jun 12 57916.2 19.35(0.05) 19.27(0.10) 19.22(0.11) 
2017 Jun 13 57917.2 19.34(0.11) 19.32(0.11) 19.29(0.14) 
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C (err) 
(mag) 


19.03(0.12) 
18.94(0.11) 
18.88(0.11) 
18.85(0.11) 
19.04(0.09) 
19.01(0.09) 
18.97(0.08) 
19.02(0.09) 


18.57(0.08) 
18.71(0.08) 
18.75(0.08) 
18.70(0.08) 


18.73(0.08) 
18.72(0.08) 
18.64(0.08) 
18.66(0.08) 
18.75(0.08) 
18.72(0.08) 
18.67(0.08) 
18.72(0.08) 


19.09(0.13) 
18.78(0.11) 
18.94(0.12) 
18.96(0.12) 
18.88(0.08) 
18.94(0.08) 
18.86(0.08) 
18.87(0.08) 
19.01(0.09) 
19.21(0.09) 
19.12(0.09) 
19.10(0.09) 


G (err) 
(mag) 


>20.7 
18.61(0.01) 
18.59(0.01) 
18.20(0.01) 
18.23(0.01) 
18.73(0.01) 
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Telescope / Instrument 


PSST 
Gaia 
Gaia 
Gaia 
Gaia 
Gaia 
Gaia 
MLS 
MLS 
MLS 
MLS 
CRTS 
CRTS 
CRTS 
CRTS 
PSST 
PSST 
PSST 
PSST 
PSST 
CRTS 
CRTS 
CRTS 
CRTS 
PSST 
PSST 
PSST 
PSST 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
PSST 
PSST 
PSST 
PSST 
MLS 
MLS 
MLS 
MLS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
CRTS 
FLWO48 
FLWO48 
FLWO48 
FLWO48 
FLWO48 
PSST 
PSST 
PSST 
FLWO48 
FLWO48 
FLWO48 
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Table 6 

(Continued) 
UT Date MJD g (err) r (err) i (err) C (err) G (err) Telescope /Instrument 
(year month day) (mag) (mag) (mag) (mag) (mag) 
2017 Jun 16 57920.2 19.53(0.12) 19.40(0.04) 19.31(0.12) MDM 
2017 Jun 18 57922.3 Ue. 19.40(0.09) 19.31(0.13) MDM 
2017 Jun 30 57934.1 19.61(0.14) Ut Un FLWO48 
2017 Jun 30 57934.2 19.54(0.09) 19.41(0.11) 19.06(0.22) FLWO48 
2017 Jun 30 57934.2 te 19.44(0.12) 19.30(0.14) FLWO48 
2017 Jul 2 57936.1 19.39(0.11) Ut Ue. FLWO48 
2017 Jul 2 57936.2 19.51(0.10) 19.44(0.11) 19.41(0.11) FLWO48 
2017 Jul 2 57936.2 Un 19.51(0.12) 19.35(0.12) Ut FLWO48 
2017 Jul 2 57936.5 Ue. Ut Ue 19.32(0.09) CRTS 
2017 Jul 2 57936.5 Ue. Ut te 19.32(0.09) CRTS 
2017 Jul 2 57936.5 te te te 19.38(0.09) CRTS 
2017 Jul 2 57936.5 te te 19.24(0.09) CRTS 
2017 Jul 08 57942.2 19.60(0.20) 19.20(0.11) Ut FLWOA4S8 
2017 Jul 08 57942.2 19.47(0.11) 19.18(0.16) Ue FLWO48 
2017 Jul 22 57956.9 te te 19.34(0.01) Gaia 
2017 Jul 22 57956.9 te te te 19.40(0.01) Gaia 
2017 Nov 15 58072.5 21.12(0.10) 20.92(0.15) 20.67(0.16) tee FLWO48 
2017 Nov 15 58072.5 20.98(0.12) 20.85(0.10) 20.45(0.15) FLWO48 
2017 Nov 19 58076.5 21.24(0.12) 20.98(0.13) 20.70(0.18) FLWO48 
2017 Nov 19 58076.5 21.25(0.12) 21.03(0.09) 20.74(0.20) FLWO48 
2017 Nov 19 58076.5 21.27(0.12) 20.95(0.15) 20.49(0.19) FLWO48 
2018 Jan 31 58149.5 Ue. 22.14(0.22) Ue. MMTCam 
2018 Jan 31 58149.5 22.23(0.10) Ut Uer FLWOA4S8 
2018 Jan 31 58149.6 Ue. Ut 21.41(0.19) MMTCam 
2018 Feb 05 58154.5 22.36(0.18) 21.77(0.27) 21.62(0.17) MMTCam 
2018 Feb 25 58174.4 22.43(0.11) 22.11(0.09) 21.50(0.10) MMTCam 
2018 Feb 26 58175.4 22.31(0.20) 22.06(0.12) Ue MMTCam 
2018 Mar 16 58193.4 EE EE 21.61(0.11) MMTCam 
2018 Mar 18 58195.2 22.38(0.25) Un te LDSS 
2018 Mar 18 58195.3 te 22.08(0.09) 21.57(0.15) LDSS 
2018 Apr 5 58213.4 22.33(0.31) 22.27(0.32) te FLWO48 
2018 Apr 5 58213.5 Ue Ut 21.64(0.32) FLWO48 
2018 May 17 58255.3 Ue 21.98(0.12) 21.56(0.11) MMTCam 
2018 Jul 09 58309.0 E Ut 21.95(0.28) LDSS 
2018 Jul 09 58309.0 Uer Ue 21.95(0.28) LDSS 
2018 Jul 09 58309.0 E Un 21.95(0.28) LDSS 
2018 Dec 09 58461.5 EE 22.73 (0.44) EE MDM 
2019 Jan 31 58514.3 Ue. 23.43 (0.08) Ue. LDSS 
2019 Jan 31 58514.3 Ue. Ut 23.52 (0.08) LDSS 


Note. Optical photometry of SN 2016iet, not corrected for Galactic extinction. Optical gri magnitudes are in the AB system, C magnitudes are calibrated to Vega 
magnitudes, and G magnitudes come from the Gaia Photometric Alerts Index. The photometric error bars are shown in parenthesis. The “>” represents a 30 upper 


limit. 
Table 7 
Optical Spectroscopy of SN 2016iet 

UT Date MJD Phase* Exposure Time Airmass Wavelength Range Telescope + Instrument Grating Slit width 

(d) (s) (A) (lines/mm) (arcsec) 
2017 Apr 30 57873.00 +132 1200 1.79 4160-9340 Baade + IMACS 300 019 
2017 May 23 57896.33 +153 1200 1.481 3310-8480 MMT + Blue Channel 300 1” 
2017 Jun 17 57921.00 +177 3600 1.70 3920-9050 MDM + OSMOS 1600 1” 
2017 Jun 29 57933.22 +187 1200 1.51 3320-8560 MMT + Blue Channel 300 1” 
2017 Nov 26 58083.51 +328 1200 1.53 3320-8590 MMT + Blue Channel 300 1" 
2018 Jan 27 58145.44 4-386 1800 1.07 3290-8550 MMT + Blue Channel 300 1" 
2018 Feb 9 58158.50 4-399 1200 1.05 3820-9200 MMT + Binospec 270 1” 
2018 Mar 17 58194.24 +433 1800 1.79 3900-10600 Clay + LDSS3c 400 1” 
2019 Mar 4 58556.58 +772 7200 1.05 5300-10200 Gemini + GMOS 400 1” 
Note. 


è Rest-frame days since peak bolometric brightness. 


19 


THE ASTROPHYSICAL JOURNAL, 881:87 (20pp), 2019 August 20 


ORCID iDs 


Sebastian Gomez © https: //orcid.org /0000-0001-6395-6702 
Edo Berger (9 https: //orcid.org /0000-0002-9392-968 1 

Matt Nicholl © https: /orcid.org/0000-0002-2555-3192 

Peter K. Blanchard O https: //orcid.org /0000-0003-0526-2248 
V. Ashley Villar © https: //orcid.org /0000-0002-5814-4061 
Joel Leja © https: //orcid.org /0000-0001-6755-1315 

Griffin Hosseinzadeh © https: //orcid.org /0000-0002- 
0832-2974 

Philip S. Cowperthwaite © https: //orcid.org /0000-0002- 
2478-6039 


References 


Aguilera-Dena, D. R., Langer, N., Moriya, T. J., & Schootemeijer, A. 2018, 
ApJ, 858, 115 

Arcavi, L, Gal-Yam, A., Sullivan, M., et al. 2014, ApJ, 793, 38 

Arcavi, L, Howell, D. A., Kasen, D., et al. 2017, Natur, 551, 210 

Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 

Astropy Collaboration, Price-Whelan, A. M., Sipócz, B. M., et al. 2018, AJ, 
156, 123 

Barbary, K. 2016, extinction, vv0.3.0, Zenodo, doi:10.5281 /zenodo.804967 

Barkat, Z., Rakavy, G., & Sack, N. 1967, PhRvL, 18, 379 

Ben-Ami, S., Gal-Yam, A., Mazzali, P. A., et al. 2014, ApJ, 785, 37 

Blanchard, P. K., Berger, E., & Fong, W.-f. 2016, ApJ, 817, 144 

Blanchard, P. K., Nicholl, M., Berger, E., et al. 2018, ApJ, 865, 9 

Boubert, D., Guillochon, J., Hawkins, K., et al. 2018, MNRAS, 479, 2789 

Branch, D., Baron, E., Thomas, R. C., et al. 2004, PASP, 116, 903 

Charlot, S., & Fall, S. M. 2000, ApJ, 539, 718 

Chatzopoulos, E., Wheeler, J. C., Vinko, J., Horvath, Z. L., & Nagy, A. 2013, 
ApJ, 773, 76 

Chevalier, R. A., & Irwin, C. M. 2011, ApJL, 729, L6 

Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 

Chomiuk, L., Chornock, R., Soderberg, A. M., et al. 2011, ApJ, 743, 114 

Chornock, R., Berger, E., Gezari, S., et al. 2014, ApJ, 780, 44 

Crowther, P. A. 2013, MNRAS, 428, 1927 

De, K., Kasliwal, M. M., Cantwell, T., et al. 2018, ApJ, 866, 72 

Dessart, L., Waldman, R., Livne, E., Hillier, D. J., & Blondin, S. 2013, 
MNRAS, 428, 3227 

Dexter, J., & Kasen, D. 2013, ApJ, 772, 30 

Drake, A. J., Djorgovski, S. G., Mahabal, A., et al. 2009, ApJ, 696, 870 

Dressler, A., Bigelow, B., Hare, T., et al. 2011, PASP, 123, 288 

Drout, M. R., Chornock, R., Soderberg, A. M., et al. 2014, ApJ, 794, 23 

Fabricant, D. G., Epps, H. W., Brown, W. L., Fata, R. G., & Mueller, M. 2003, 
Proc. SPIE, 4841, 1134 

Foreman-Mackey, D. 2016, JOSS, 1, 24 

Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 
125, 306 

Fransson, C., & Chevalier, R. A. 1989, ApJ, 343, 323 

Galama, T. J., Vreeswijk, P. M., van Paradijs, J., et al. 1998, Natur, 395, 670 

Gal-Yam, A., Mazzali, P., Ofek, E. O., et al. 2009, Natur, 462, 624 

Gelman, A., & Rubin, D. B. 1992, StaSc, 7, 457 

Gezari, S., Heckman, T., Cenko, S. B., et al. 2009, ApJ, 698, 1367 

Guillochon, J., Nicholl, M., Villar, V. A., et al. 2018, ApJS, 236, 6 

Guillochon, J., Parrent, J., Kelley, L. Z., & Margutti, R. 2017, ApJ, 835, 64 

Heger, A., Woosley, S., Baraffe, L, & Abel, T. 2002, in Proc. MPA/ESO/ 
MPE/USM Joint Astronomy Conf., Lighthouses of the Universe: The Most 
Luminous Celestial Objects and Their Use for Cosmology, ed. M. Gilfanov, 
R. Sunyeav, & E. Churazov (Berlin: Springer-Verlag), 369 

Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 

Hony, S., Gouliermis, D. A., Galliano, F., et al. 2015, MNRAS, 448, 1847 

Hook, I. M., Jørgensen, I., Allington-Smith, J. R., et al. 2004, PASP, 116, 425 

Huber, M., Chambers, K. C., Flewelling, H., et al. 2015, ATel, 7153, 1 

Hunter, J. D. 2007, CSE, 9, 90 


20 


Gomez et al. 


Inserra, C., Smartt, S. J., Jerkstrand, A., et al. 2013, ApJ, 770, 128 

Jerkstrand, A. 2017, in Handbook of Supernovae, ed. A. W. Alsabti & 
P. Murdin (New York: Springer), 795 

Jerkstrand, A., Smartt, S. J., & Heger, A. 2016, MNRAS, 455, 3207 

Karamehmetoglu, E., Taddia, F., Sollerman, J., et al. 2017, A&A, 602, A93 

Kasen, D., & Bildsten, L. 2010, ApJ, 717, 245 

Kasen, D., Woosley, S. E., & Heger, A. 2011, ApJ, 734, 102 

Kasliwal, M. M., Kulkarni, S. R., Gal-Yam, A., et al. 2012, ApJ, 755, 161 

Kelly, P. L., & Kirshner, R. P. 2012, ApJ, 759, 107 

Kennicutt, R. C. J. 1998, ARA&A, 36, 189 

Kobulnicky, H. A., Kennicutt, R. C. J., & Pizagno, J. L. 1999, ApJ, 514, 544 

Leja, J., Johnson, B. D., Conroy, C., van Dokkum, P. G., & Byler, N. 2017, 
ApJ, 837, 170 

Leloudas, G., Chatzopoulos, E., Dilday, B., et al. 2012, A&A, 541, A129 

Lunnan, R., Chornock, R., Berger, E., et al. 2014, ApJ, 787, 138 

Lunnan, R., Chornock, R., Berger, E., et al. 2015, ApJ, 804, 90 

Lunnan, R., Chornock, R., Berger, E., et al. 2016, ApJ, 831, 144 

Lunnan, R., Fransson, C., Vreeswijk, P. M., et al. 2018, NatAs, 2, 887 

Lunnan, R., Kasliwal, M. M., Cao, Y., et al. 2017, ApJ, 836, 60 

Martini, P., Stoll, R., Derwent, M. A., et al. 2011, PASP, 123, 187 

Mazzali, P. A., Moriya, T. J., Tanaka, M., & Woosley, S. E. 2019, MNRAS, 
484, 3451 

Milisavljevic, D., Patnaude, D. J., Raymond, J. C., et al. 2017, ApJ, 846, 50 

Modjaz, M., Bianco, F. B., Siwek, M., et al. 2019, arXiv:1901.00872 

Modjaz, M., Liu, Y. Q., Bianco, F. B., & Graur, O. 2016, ApJ, 832, 108 

Montero-Dorta, A. D., & Prada, F. 2009, MNRAS, 399, 1106 

Moriya, T. J., Nicholl, M., & Guillochon, J. 2018, ApJ, 867, 113 

Nadyozhin, D. K. 1994, ApJS, 92, 527 

Nicholl, M. 2018, RNAAS, 2, 230 

Nicholl, M., Berger, E., Blanchard, P. K., Gomez, S., & Chornock, R. 2019, 

ApJ, 871, 102 

Nicholl, M., Berger, E., Margutti, R., et al. 2016a, ApJL, 828, L18 

Nicholl, M., Berger, E., Smartt, S. J., et al. 2016b, ApJ, 826, 39 

Nicholl, M., Blanchard, P. K., Berger, E., et al. 2018, ApJL, 866, L24 

Nicholl, M., Guillochon, J., & Berger, E. 2017, ApJ, 850, 55 

Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2013, Natur, 502, 346 

Nicholl, M., Smartt, S. J., Jerkstrand, A., et al. 2015, ApJL, 807, L18 

Oliphant, T. E. 2007, CSE, 9, 10 

Owocki, S. P. 2015, ASSL, 412, 113 

Perets, H. B., Gal-Yam, A., Mazzali, P. A., et al. 2010, Natur, 465, 322 

Quimby, R. M., De Cia, A., Gal-Yam, A., et al. 2018, ApJ, 855, 2 

Quimby, R. M., Kulkarni, S. R., Kasliwal, M. M., et al. 2011, Natur, 474, 487 

Rakavy, G., & Shaviv, G. 1967, ApJ, 148, 803 

Rezaeikh, S., Javadi, A., Khosroshahi, H., & van Loon, J. T. 2014, MNRAS, 
445, 2214 

Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 

Schmidt, G. D., Weymann, R. J., & Foltz, C. B. 1989, PASP, 101, 713 

Science Software Branch at STScI 2012, PyRAF: Python alternative for IRAF, 
Astrophysics Source Code Library, ascl:1207.011 

Smithsonian Astrophysical Observatory 2000, SAOImage DS9: A utility for 
displaying astronomical images in the X11 window environment, 
Astrophysics Source Code Library, ascl:0003.002 

Stevenson, K. B., Bean, J. L., Seifahrt, A., et al. 2016, ApJ, 817, 141 

Stritzinger, M., Mazzali, P., Phillips, M. M., et al. 2009, ApJ, 696, 713 

Taddia, F., Sollerman, J., Fremling, C., et al. 2019, A&A, 621, A64 

Valenti, S., Elias-Rosa, N., Taubenberger, S., et al. 2008, ApJL, 673, L155 

van der Walt, S., Colbert, S. C., & Varoquaux, G. 2011, CSE, 13, 22 

Vreeswijk, P. M., Leloudas, G., Gal-Yam, A., et al. 2017, ApJ, 835, 58 

Woosley, S. E. 2010, ApJL, 719, L204 

Woosley, S. E. 2017, ApJ, 836, 244 

Woosley, S. E., Blinnikov, S., & Heger, A. 2007, Natur, 450, 390 

Wyrzykowski, L., Hodgkin, S., Blogorodnova, N., Koposov, S., & Burgon, R. 
2012, in Proc. GAIA-FUN-SSO 2012 Conf., 2nd Gaia Follow-up Network 
for Solar System Objects, ed. P. Tanga & W. Thuillot (Paris: Observatoire 
de Paris), 21 

Yan, L., Lunnan, R., Perley, D. A., et al. 2017, ApJ, 848, 6 

Yoshida, T., Umeda, H., Maeda, K., & Ishii, T. 2016, MNRAS, 457, 351 

Yusof, N., Hirschi, R., Meynet, G., et al. 2013, MNRAS, 433, 1114 


